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We consider 2D Dirac fermions in the presence of three types of disorder: random scalar potential, 
random gauge potential and random mass with long-range correlations decaying as a power law. 
Using various methods such as the self-consistent Born approximation (SCBA), renormalization 
group (RG), the matrix Green function formalism and bosonisation we calculate the density of 
states and study the full counting statistics of fermionic transport at lower energy. The SCBA 
and RG show that the random correlated scalar potentials generate an algebraically small energy 
scale below which the density of states saturates to a constant value. For correlated random gauge 
potential, RG and bosonisation calculations provide consistent behavior of the density of states 
which diverges at zero energy in an integrable way. In the case of correlated random mass disorder 
the RG flow has a nontrivial infrared stable fixed point leading to a universal power-law behavior of 
the density of states and also to universal transport properties. In contrast to uncorrelated case the 
correlated scalar potential and random mass disorders give rise to deviation from the pseudodiffusive 
transport already to lowest order in disorder strength. 

PACS numbers; 73.63.-b,73.22.Pr,73.23.-b 



I. INTRODUCTION 

In the recent years significant attention has been 
attracted by materials exhibiting two dimensional 
fermionic excitations with linear dispersion relation near 
the Fermi level. These excitations share many proper- 
ties with massless relativistic particles but with a ve- 
locity reduced with respect to the speed of light. The 
seminal example of such material is graphenCfi the low 
energy properties of which are described by two dimen- 
sional (2D) gas of Dirac fermions More recently, 2D 
Dirac fermions have also emerged as the effective low 
energy degree of freedom in the surface states of 3D 
topological insulators fi"— such as the materials in the 
Bi2Se3 family^ and strained HgTe.* Dirac excitations 
have been also found in unconventional superconductor 
with d-wave symmetry^"^^, in the quasi-2D organic con- 
ductor a-(BEDT-TTF)2 I3 under pressurcr^"— and even 
in photonic crystals. 

The peculiar features of Dirac fermions lead to unfa- 
miliar transport properties in all these materials. Of par- 
ticular interest are the vanishing of the density of states 
at the Dirac (or neutrality) point, and the unusual scat- 
tering properties of the Dirac particles, a striking exam- 
ple of which being provided by the so-called Klein tun- 
neling phenomenon: an excitation incident normal to a 
potential barrier crosses this barrier completely even for 
energies smaller than the barrier height,- A consequence 
of that is the total absence of backscattering for these 
Dirac fermions. The peculiarities of the spectrum near 
the Dirac point also affect transport strongly. In undoped 
graphene, due to the evanescent nature of the states at 
the Dirac point, transport in a clean sample is similar 
to that in a diffusive wire.^" Moreover, all the different 
cumulants of current fluctuations in a graphene sample 
behave as in a diffusive wire rather than in a ideal metal- 
lic system. In particular, pseudodiffusive conductance 



scales with the length of the system L as 1/L in contrast 
with the L° scaling of ideal conducting systems. The 
associated Fano factor defined as the ratio between the 
shot noise power and the current has the same universal 
value F = 1/3 as for diffusive metallic wires instead of 
F = for ideal conductors. Remarkably, the conductiv- 
ity minimum of graphene at the neutrality point due to 
evanescent modes is of order e^/h, i.e. finite despite the 
vanishing of the density of states. Though the conduc- 
tivity minimum remains almost constant in very broad 
temperature range its sample-dependance indicates the 
importance of disorder for the transport properties of 
graphene 1^ 

Due to these specificities, numerous studies have focus 
on the effect of a random scattering potential on trans- 
port of Dirac states.—"— Indeed, various kinds of dis- 
order are naturally present in real materials, and affect 
in a dominant way the electronic transport properties. 
They can be of different origin: lattice defects, impu- 
rities, ripples in graphene sheet that distort locally the 
lattice, adatoms deposited on the surface of a graphene 
sampler^ atomic steps on the surface of topological insu- 
lators, etc. Theoretical research on disordered 2D Dirac 
fermions was also motivated initially by its relevance to 
Quantum Hall transitions*^ Building on this pioneering 
work, recent studies have focused on the effect of the dif- 
ferent types of disorder on the transport properties of 
Dirac fermions as in the limit of low energy, i.e. around 
the Dirac point, as well as away from half filling. In 
the present paper, we do not consider the highly doped 
weak localization regime^^S? corresponding to kplo ^ 1, 
where Iq is the mean free path, and concentrate mostly on 
the transport near the Dirac point where kp ~> 0. It has 
been shown that the conductivity at half filling depends 
not only on the type of disorder but also on the infrared 
cutoff, so that it potentially can depend on the geometry 
of the physical setup fSi The role of infrared cutoff can 
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be played by either the mean free path, the Fermi length 
or the size of the system. These different cases allow 
to identify several transport regimes Fixing geometry, 
for example to wide-and-short rectangle with many prop- 
agating transverse modes one can compute the conduc- 
tance and the Fano factor for each regime;^ Most of the 
previous studies addressed the case where the scattering 
potential is uncorrelated in space. This is the case for 
instance if it is originated from localized point-like scat- 
terers distributed independently from each other. How- 
ever, several physically relevant types of disorder sources 
exhibit long-range correlations. 

For example, a graphene sheet is known to develop 
static shape fluctuations due to the unavoidable thermo- 
dynamic instability of 2D crystals with respect to both 
crumpling and bending. These ripples survive at low 
temperatures and can be viewed as a static random gauge 
potential playing the role of a quenched disorder on the 
electronic time scale (see below) . The theory of 2D elas- 
tic membranes predicts the strength of the local height 
fluctuations which give rise to long-range algebraic cor- 
relation of this random gauge potentia h^^i^"" — A second 
example, in the case of the surface states of topologi- 
cal insulators, is surface roughness-^ which is one of the 
dominant form of disorder in these materials. A typical 
roughness created by atomic steps can induce a scattering 
potential with algebraically decaying correlations.-^''^^ As 
a last example, the adsorption of magnetic adatoms on 
the surface of topological insulators has been proposed 
as a way to control the electronic properties of the sur- 
face states.— If the characteristic spin flip time of mag- 
netic adatoms exceeds the mean free time of the elec- 
trons in the surface states these adatoms can be also 
viewed as a source of quenched disorder of both random 
gauge and random mass types. In the vicinity of the 
paramagnetic-ferromagnetic transition, induced by the 
RKKY-type surface interactions, critical magnetization 
fluctuations will give rise to a quenched disorder on elec- 
tronic time scales with power law correlations in space. 

Motivated by these physical examples, in this paper we 
consider the general properties of 2D Dirac fermions in 
the presence of various weak random potentials possess- 
ing algebraic spatial correlations. Using several analyti- 
cal techniques, we consider perturbatively these disorder 
potentials. We focus on the effect of these long-range cor- 
relations on the density of states and also transport prop- 
erties discussing the cases relevant for the three above 
mentioned examples. In particular, we will compute the 
unknown to our knowledge density of states for corre- 
lated random potential and random mass. We will show 
that the previous estimation^ of the density of states for 
correlated random gauge potential is wrong. We will find 
the correct density of states using two different methods: 
renormalization group and bosonisation. We will develop 
a framework to study transport properties in the pres- 
ence of correlated disorder using the matrix Green func- 
tion formalism introduced by Nazarovj^^i^ Our approach 
goes beyond the previous work of Khveshchenko ^^'^^ who 



also considered LR correlated potentials but focused on 
the multi fractal spectrum of wave functions at the Dirac 
points and the conductance within the self consistant 
Born approximation (SCBA). 

The paper is organized as follows. Section |lT] intro- 
duces the model. In Sec. IIIII we consider the SCBA ap- 
proximation. In Sec. IIVI using the matrix Green func- 
tion formalism we study the full counting statistics for a 
wide-and-short rectangle sample at the neutrality point. 
In Sec.|V]we derive RG equations to one- loop order and 
discuss the properties of the systems with different type 
of disorder. In Sec. I VII we use bosonisation technique for 
systems with LR correlated random gauge potential. In 
Sec. IVIII wc summarize the obtained results. 

II. MODEL 

A. Single flavor Dirac model 

Whenever they appear in a 2D or quasi-2D material, 
like graphene, lattice Dirac fermions are constrained by 
the Nielsen-Ninomiya theorem^ to appear by pair of 
species. Practically, this implies the existence of an even 
number of Dirac cones in the first Brillouin zone when 
considering the low energy dispersion relation. Indeed, in 
graphene two Dirac cones exist at the inequivalent points 
K and K' = —K at the zone boundary. However, any 
potential varying on scales much larger that the atomic 
scale Ti/K will leave the two Dirac cones uncoupled. In 
this case, the effect of the potential can be described by 
considering its effect on a single Dirac cone, treating the 
presence of the other Dirac point as an effective degen- 
eracy. We thus lead to consider a single species of non- 
interacting massless 2D Dirac fermions in the presence of 
a random potential, described by the Hamiltonian 

H = Ho + V{x,y), (1) 

where Hq is the kinetic Hamiltonian of free Dirac 
fermions with the Fermi velocity vq, 

Ho = -ivoia^dr, + Uydy), (2) 

and V{x, y) is a random disorder potential. Here and in 
the following ctq — 1, cr^, /J. — x,y,z are the respective 
Pauli matrices, and we set h = 1 for convenience. This 
type of potential, without any Fourier component which 
couple the different Dirac species, is often denoted as a 
long-range potential. This notation, which refers to cor- 
relations at the scale of the lattice space, should not be 
confused with the long-range (LR) correlation in space 
on which we focus in this paper. The latter characterizes 
the long distance g ~ behavior of the random poten- 
tial correlations. Note that in the case of graphene, as 
well as in the quasi-2D a-(BEDT-TTF)2 I3, the Pauh 
matrices entering the relativistic kinetic Hamiltonian re- 
fer to a pseudo-spin describing the relative weight of the 
electronic wave function on two sub-lattices. Thus, the 
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coupling of the random potential to these two sub lat- 
tices will be reflected in the parametrization introduced 
below for this potential in terms of these Pauli matrices. 

The case of surface states of topological insulators is 
differential^ In these materials, a strong spin orbit in- 
teraction opens a gap for bulk states. The non trivial 
topological order characterizing the filled bands of this 
insulator implies the existence of Dirac fermions surface 
states. Since they are not constrained by the Nielsen- 
Ninomiya theorem, they occur around an odd number 
of Dirac points in the first Brillouin zone. In the sim- 
plest topological insulators, a single Dirac cone exists at 
the surface of these insulators, the properties of which 
control the surface transport properties of the material. 
In this case, the relativistic kinetic Hamiltonian of the 
surface states contains a real magnetic term reflecting 
the bulk spin orbit interaction. This term should be an 
odd function of the electron spin, which we take for sim- 
plicity as the in plane component of the spin, obtaining 
effectively Eq. (HI). In this case, the parametrization of 
disorder depends on the coupling of the corresponding 
potential to the spin of the electrons. In topological in- 
sulators, the bulk topological order at the origin of this 
odd number of Dirac species, also prevents time-reversal 
invariant disorder from localizing these states. This ro- 
bustness property is in fact the result of an odd number 
of Dirac species as opposed to an even number, and will 
play no role in the decoupled cones treatment of disorder 
that we will perform. 

As discussed above, we parametrize the disorder po- 
tential using the following decomposition 



(3) 



In the usual terminology of disordered graphene the term 
with /i = is called random potential disorder, the terms 
with fj, = x,y random gauge disorder and the term with 
fi = z random mass disorder. We will keep this terminol- 
ogy, even though for systems with real spin (instead of 
pseudo-spin) like topological insulators they may have a 
different physical interpretation. Fox instance, in topo- 
logical insulators the terms with fj, = x,y, z in Eq. ^ 
correspond to real random magnetic impurities on the 
conducting surface. 

In what follows the disorder potentials V"^(r) [r ~ 
{x, y}] are taken to be random and Gaussian with 
(V^(r)) = and correlators 



(y4r)K(r')) = 27ri;oV.9M(r - r'). 



(4) 



For the sake of convenience we fix the form of the corre- 
lator in Fourier space 

(V;(k)K(k')) - {2TT)\iS{k + k')S^,{a^ +f3^\kr-^).{5) 
This form of the correlator corresponds in real space to 



The exponent a is determined by the nature of disorder 
correlations or by internal or fractal dimension of the ex- 
tended defects. In the presence of extended defects of in- 
ternal dimension Sd randomly orientated the correspond- 
ing exponent a = 2 — ed- For instance, the presence linear 
dislocations or atomic steps with random orientations on 
the surface of topological insulator (e^ = 1) leads to LR 
correlated disorder with a = 1.^^ In the case of ripples of 
a graphene sheet, an evaluation of the correlation of the 
generated gauge potential provides an exponen t^^i^"" — 
a « 1.6. In the case of magnetic adatoms deposited at 
the surface of topological insulators,''^ the induced mag- 
netic disorder correlations will decay as a power law with 
a — rj where rj = 1/4 is the critical exponent describ- 
ing the magnetization correlation function in 2D Ising 
system. Hence, these examples provide strong motiva- 
tion to consider the effect of these algebraic correlations 
beyond the standard case with short-range correlation, 
like (5-corelation formally corresponding to a > 2 in two 
dimensional systems. 



B. The effective replicated action 



For non-interacting fermions, we can write the parti- 
tion function asi^ 



Z = / Dip{r,T)D'tp{r,T)e' 



(7) 



with the action given by 



S[i:,^]^ dr j d?vi^{v,T){dr-li + H)i^{v,T), (8) 

where ip and ip are anticommuting Grassman variables 
satisfying the antisymmetric boundary condition '0(r, r-t- 
/?) = — '0(r, t) and V'(r, r + /?) = — '0(r, r), and ^ is the 
chemical potential. For a time independent Hamiltonian, 
we can introduce the Fourier series decomposition: 



(9) 



where f„ — (2n -I- 1)it/ (3. This allows one to factorize 
the partition function in such a way that all terms with 
different values of are decoupled: 



Z = ]^Z(w„), 

n 



Det(zt'„ + fi — H), 
1 



(10) 



(11) 



Hence, for each value of Vn we have a two-dimensional 
time independent field theory in which the Matsubara 
frequency plays the role of a mass term. Taking the 
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zero temperature limit /3 — >■ cx3 one converts the sums 
over n into integrals over v so that we end up with 



Z = / VipVip exp 



(13) 



We are interested mostly in the properties of the undoped 
system with the Fermi energy near the Dirac cone so that 
we put /i = in what follows. Using the replica trick^ 
we derive the replicated action for the fermions at a given 
energy e = iv. To that end we introduce n replicas of the 
original system and average their joint partition function 
over disorder we obtain the effective action 



Sis) = E / 

a=l •' 



X (r)cr^i/'Q (r)] [i>fs (r')o-A''/'/3 (r')] • 



(14) 



The properties of the original system with quenched dis- 
order are then obtain by taking the limit n ^ Q. 



III. SELF-CONSISTENT BORN 
APPROXIMATION 



Let us first consider the self-consistent Born approxi- 
mation (SCBA) which is applicable only in the limit of 
weak scattering. The SCBA has been widely used to 
study the effect of uncorrelated disorderi^M^dl as well as 
the effect of Coulomb impurities^ on the Dirac fermions. 
The retarded and advanced Green functions can be ex- 
pressed in terms of the self-energy via the Dyson equation 



G(e) = Go(£)4- Go(e)S(£)G(£), 



(15) 



where Go(£, k) = (e — wock)^^ is the bare Green function 
and the dressed Green function can be written in terms 
of the self-energy as follows 



G(£,k) 



e — S(£, k) + UQcrk 



(16) 



The Green function averaged over disorder within the 
SCBA is shown schematically in Fig. [TJ One takes into 
account only the one-loop diagram contributing to the 
self-energy with the bare Green function replaced by the 
dressed one. The self consistent equation for this self en- 
ergy is simplified by its independence on the external mo- 
menta. It is easy to see that one can treat different types 
of SR and LR correlated disorder on the same footing by 
introducing the effective couplings a = aQ-\-ax-\-ay + az 
and j3 = /3o + Px + Py + Pz- The one- loop self-energy 
diagram with the dressed Green function is given by the 
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FIG. 1: The averaged over disorder Green function in the 
SCBA. 



integral 



S(e)= / 27r4g(fc)Go(£,k) 



A/-UO 



A(e) 



v1g{k)kdk 
X^{e)~vik^'' 



(17) 



where we have introduced the UV momentum cutoff 
A/vq. The function X{e) = e — S(£) has to be deter- 
mined self-consistently. Once Eq. PT)) is solved the den- 
sity of states can be computed using the retarded Green 
function as follows 



Pis) 



-Im Tr / G^(e,k) 



ImA:(£)ln 



A2 



XHe) 



(18) 



We now consider separately the cases of the SR and LR 
correlated disorder. 

SR correlated disorder. In this case the disorder cor- 
relator reads from eq. ^ g{k) = a so that the SCBA 
equation (fT7|) reduces to 



X{e) = e+-X{e)\n 



X^ie) 



(19) 



The solution of Eq. ([TS]) has two branches which cor- 
respond to the retarded and advanced Green functions. 
They can be written explicitly in terms of the Lambert 
function'*^ W{x) as follows^ 



X{e)^e/{aWi±ie/iaTo))). 



(20) 



Here and below, the upper sign corresponds to retarded 
and the lower sign to advanced functions. The weak dis- 
order introduces a new exponentially small energy scale 
Fo = Ae^i/". For £ < Fq one can expand Eq.lgO]) in 
small e, this yields X{0) = ^iTo + e/a. At the Dirac 
point the self-energies are pure imaginary S(0) — T*ro- 
This would imply a finite density of states at the Dirac 
cone that contradicts with the results obtained by dif- 
ferent methods. For instance, for uncorrelated random 
gauge disorder one expects p{e) — e^l^^^ with z — \ + a. 
in weak disorder case^^ (a < 2) and z — (8a) — 1 in 
strong disorder case**^ [a > 2). It was argued in Ref. [s^ 
that the failure of SCBA in the vicinity of the Dirac cone 
for e < Fq is due to importance of the diagrams with 
crossed disorder lines, neglected within the SCBA which 
takes into account only the non-crossed ones. Another 
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reason for the failure of the SCBA is the divergence of 
the fermion wavelengths at the Dirac point rendering the 
SCBA uncontrollable, i.e. without a small parameter. 
Nevertheless one can still rely on the SCBA for energies 
e 3> Fq. For Fq ^ e ^ A and a <^ 1 one can derive an 
approximate solution of ([T9l) by iteration of the Eq. ([T9| . 
To lowest order one obtains2^ 









A" 


^1 + aln— ^ 


± —TTae 


1 - 


h 2a In — 




2 




s _ 



(21) 



Reexpressing X \n{-A'^ / X^) in Eq. ^ using Eq. (jl9p 
we get the density of states 



Pscba(£) 



-Im ■ 



1 



n^v^a^ "' W{e/{iaro)))' 
For e ^ Fq this simplifies to 



Pscba(£) 



A' 

1 + 2aln — 

e 



(22) 



(23) 



while below the energy scale Fq the density of states sat- 
urates at a finite value which is in fact correct only for 
random potential disorder. ^'^ 

LR correlated disorder. The disorder correlator having 
the form g{k) = l3k''~^ (see eq. jS])) yields the SCBA 
equation of the following form 



X{s) 



a^,2-a 




aX{e) ' ' V'2' ' 2'X2(£) 



, (24) 



where 2Fi{a,b,c;z) is the Gauss hypergeometric 
function.— In the limit a — ?► 2 we recover Eq. ((T9)) with 
a replaced by p. Similarly to the previous case Eq. ([M)) 
has two solutions corresponding to the retarded and ad- 
vanced Green functions. At the Dirac point the self- 
energies are pure imaginary: E(0) = T*F2-a and disorder 
induces a generalized small energy scale 



F2-a = A 



Ae-1/^ 



'2sin(^) {2~a + P) 



{2-a)Trl3 
1 



(2 - a) ^ + — 




0((2-a)2) 



,(25) 



where we have introduced the dimensionless disorder 
strength /3 — /3(A/wo)"~^ and in the second line we 
have performed an expansion to the first order in 2 — a. 
For finite 2 — a the energy scale F2-a is only alge- 
braically small in disorder at variance with the expo- 
nentially small energy scale Fq for uncorrelated disorder. 
Iterating Eq. (IMl) one obtains an approximate solution, 
valid for F2_a ^ e ^ A and /? ^ 1. To lowest order 
that solution reads 



X{e) = s{l + PU{£)) ± i^-^h [1 + 2^C/(£)] , (26) 



, , vr / A\ ^ " 7ra 1 

^(^)--2 7 ^°*T-23^' 



(27) 



In the limit of a — )■ 2 one obtains U{£) ~ ln(A/£). Substi- 
tuting the solution ([25]) in Eq. (ITS)) we obtain the density 
of states for F2-a ^ e 



Pscba(£) 



2t:vI 



1 + ;3 ( In - + [/(e) 



(28) 



which reproduces the density of states ([25)) in the limit 
of SR correlated disorder a ^ 2. 



IV. FULL COUNTING STATISTICS 

A. Matrix Green function formalism 

We now consider the transport properties of 2D Dirac 
fermions propagating in a rectangular sample of size 
L X W . The schematic setup is shown in Fig. [5| We 
assume that the perfect metallic leads are attached to 
the two sides of the width W with the distance L <^ W 
between them. We model the leads as heavily doped re- 
gions described by the same Dirac Hamiltonian but with 
the chemical potential e_F 3> e shifted far from the chem- 
ical potential e ~ in the bulk which is close to the 
Dirac point. The large number of propagating modes in 
the leads are labeled by the momentum p„ = 2Tm/W 
in y direction with n = 0, ±1,.. ± epW/{2TTVQ). In the 
limit W ^ L in which many modes A'^ ^ 1 contribute 
to transport one can neglect the boundary conditions 
at y = ±W/2 and treat p„ as a continuous variable 
p. It is convenient to switch from the coordinate rep- 
resentation to the mixed channel-coordinate representa- 
tion ip{x,y) — ^ ipnix) where n enumerates the transverse 
modes. Using this basis one can describe the wave func- 
tions in the leads by two vectors c'" = [{a+},{6~}] and 
c°"* = [{a~}, {&^}] where a„ and 6„ refer to the ampli- 
tude of waves in the left and in the right lead, respec- 
tively. The sign " refers to the waves moving to right 
and the sign "-" to the waves moving to left. These 
two vectors are related by the scattering matrix S as 
c°"* = 5c"\ In the lead subspace it has the standard 
structure)^ 



f t' 
i f' 



(29) 



where we use the "hat" notation for matrices defined in 
the channel space. The conservation of particles implies 
that iS is a unitary matrix and that the four Hermitian 
matrices ii^ , i'H'^ , 1 — ff^ , 1 — f'^f'^ have the same set of 
eigenvalues r„ each of them is a real number between 
and 1. The transport statistics is completely determined 
by the matrix of transmission amplitudes tmn between 
channels m and n in the leads since the transmission 
probabilities of the system are given by the eigenvalues 
T„ of the matrix iH^ 

The scattering matrix S relates incoming to outgoing 
states. An alternative formulation is based on the the 
transfer matrix T which relates the states in the left lead 
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FIG. 2; (Color online) The setup for two-terminal transport 
measurements: a 2D disordered sample of size W x L with 
perfect leads attached on opposite sides. 



to states in the right lead: c"^^^ = Tc'°^*. The waves in 
the left and the right leads are given by the vectors c'*^^* = 
[{a+},{a-}] and c^s^t = [{6+}, {6"}], respectively. One 
can show that the eigenvalues of TT^ appear in pairs of 
the form e^^^" with A„ > related to the transmission 
eigenvalues by T„ = 1 / cosh^ A„ . 

The transmission eigenvalues allow one to calculate 
a variety of transport properties. In the limit of large 
number of channels one can introduce the distribution 
function P{T). By definition / dTP{T) gives the total 
number of open channels. The first two moments of the 
distribution give the Landauer conductance 



p2 p2 foo 

G=—TiiH^— dTTPm. (30) 
h h ' 







and the Fano factor 

Tr(£t£)2 
Trdt 



F = 1 



= 1 - 



dTT^P{T)^ j [J^ dTTP{T)^ , (31) 



which describes the power spectrum of the noise due to 
discreetness of the charge carriers at zero frequency and 
average current /: Pq = 2eFI. Note, that for graphene 
one has to multiply Eq. (|30|) by the factor of 4 accounting 
for the spin and valley degeneracy. It is also convenient 
introduce the probability density 7^ (A) of the parameter A 
defined by T = 1/ cosh^(A), which is naturally completely 
equivalent to P{T). 

In general one can write down an integral equation for 
the transfer matrix with a kernel which depends on a 
particular realization of disorder. Iterating the integral 
equation and averaging over disorder one can compute 
the transfer matrix as an expansion in small disorder. 
However, in the case of LR correlated disorder the forth- 
coming problem of computing the transmission eigenval- 
ues seems to be a formidable task. Fortunately, there is 
an alternative way which allows one to relate P{T) to 
the free energy of an auxiliary field theory. This method 
is based on the matrix Green function formalism intro- 
duced by Nazarovi^ Instead of P{T) the statistics of the 
transmission eigenvalues can be encoded in the generat- 



ing function 

oo 

J^{z) = J2 2""^Tr(F£)" = Tr[t-i?-i - z]-\ (32) 

n=l 

All moments of P{T) can be computed using the series 
expansion of J-{z) at z = 0. The function J-{z) is regular 
in vicinity of z = and has a brunch cut along the real 
axis going from 1 to oo. Both functions are related by 
the Ricmann-Hilbert equation 



^ P{T)dT 



(33) 



and its solution is given by the jump of T{z) across the 
brunch cut 



P{T) 



1 



27riT2 



[F{l/T + - F{l/T -iQ)]. (34) 



To calculate T{z) we now adopt the matrix Green func- 
tions approach originally developed in Ref. l39l and ap- 
plied to Dirac fermions in graphene with uncorrelated 
disorder in Ref. The coefficients of the series expan- 
sion of the generating function at z = can be expressed 
in terms of the Green functions of the system as 

Tr(£t{)" ^ Tr[^G^(x,x')C'G«(x',x)]^^o^,,^^. (35) 

Here v,j. = CTj,! is the velocity operator and the retarded 
and advanced Green functions in the channel-coordinate 
representation read 

{e~H± i{))G^-^{x, x') = 5{x ~ x^)i. (36) 

The generating function can be written as a trace of an 
auxiliary two component Green function defined in the 
retarded-advanced (RA) space in the presence of ficti- 
tious field zS^ The matrix Green function is given by 



(37) 



where we use "check" notation for the objects defined in 
RA space and the operator K reads 



k(x) = 



e- H + iO -y/zvS{x) 
-^v5{x-L) e-H-iO 



(38) 



Considering the field z as a small perturbation we can 
rewrite Eq. (j37p in an integral form as follows 



G{x,x') — Gq{x,x') 

/ dxiGo{x,xi)VxiG{xi,x') (39) 



(40) 



with the kernel and the inhomogeneity given by 

/ v6{x) 
V vS{L - x) 
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One can then compute the generating function using 
1 



F{z) 



2^1 



dx Tr[V^G(x, x)] 



(41) 



that can be checked by iterating Eq. p9|) and substituting 
in Eq. (HI]) . One can relate T{z) to the object which 
plays the role of the free energy in the corresponding 
field theory. Let us rewrite the Green function (|37)) in 
coordinate representation using a functional integral over 
Grassmann variables "0 and if) 



G(r,r') = Vi!Vi!ij{Y)i;{Y')e'^ 
with the bilinear in ifj and ^ action 

I (fT[^j{v)k^{v)]. 



(42) 



(43) 



The corresponding partition function and the free energy 
can be written as follows 



Z{z) = VieiK 



VipVip e~ 



n{z) = liiZ. (44) 



It is convenient to rewrite the free energy in terms of the 
angle (f> defined by z = sin^(0/2). Direct inspection of 
Eq. shows that 



Hz) 



dVL{z) 
dz 



2 dn{(j)) 



sin (p dcf) 



(45) 



-2 arcsin ^/z 



The distribution of transmission eigenvalues V{\) can be 
calculated using the following relation 



(46) 



4>=-iT+2i\ 



from which one can easily derive P(T). Therefore, we 
have four equivalent descriptions of the transport prop- 
erties in terms of one of the following functions: P{T), 
■p(A), T{z) or fl{(j>). Any of these functions can be used 
for computing the conductance or the Fano factor. For 
instance, using the free energy one can derive the expres- 
sions 



2e2 

G=^f^"(0), 



^ _ 1 _ 2 »(iv)(0) 
^ 3 ^ 3 n"io) ' 



(47) 



where the derivatives are taken at (p = 0. Note that in 
the case of graphene the conductance ([47]) is given as 
expected per Dirac species, i.e. per spin and per valley. 



B. Expansion in disorder and diagrammatics 

In what follows we restrict our consideration to trans- 
port around the Dirac cone e = 0. The action (^5)) for 
the system including the metallic leads can be calculated 




(b) 



(c) 



(d) 



(e) 



FIG. 3: The one-loop diagrams contributing to (a) free en- 
ergy; (b) propagator; (c)-(e) disorder renormalization. 



using the kernel ([55]) with Hamiltonian ^ in which the 
free part is modified to 



Ho = -fi{x) - ivoicTjcdx + (yydy). 



(48) 



Here ^{x) — for < x < L and -l-oo otherwise accounts 
for the leads with very high chemical potential. Above we 
have treated the auxiliary field z as a perturbation. Here 
we split the kernel ((38)) into the free part including the 
auxiliary field and the interaction part: K — Kq + Ky, 
where Kq is computed using Eq. (j48|) and Ky is diagonal 
in RA space. 

We are now in the position to average the free energy 
over disorder. To that end we use the replica trick and 
introduce n copies of the original system. Performing 
averaging over disorder we obtain the replicated action 
in the following form 

n ^ 

n 3 

E E/^'r/^'^V(r-r') 



a,fi=l ti=0 

x[V;o(r)E^Va(r)][V;^(r)S^V'/3(r)] 



(49) 



where Kq is given by Eq. (|38l) with Hamiltonian (|48| and 
we defined the matrix = 1 cr^. In what follows we 
set vq = 1 unless it is written explicitly. The bare Green 
function corresponding to Kq in Eq. ([37| can be written 
in coordinate representation as 



Go{x,x';y) = 



1 



4Lcos((/)/2) 



COS 0^- 



2 Sin Tra;^ t sin -ttXq 

cos <f){^—xl) cos4>{^—x^) 



sin 4>{1 — Xq) 

zrr- 



l Sm TTX^ 

c 1 n rA'V* ^ 



I Sin 7rx -[ 
sin (px^ 
sin -Tra^Q 
sin (l)X^ 

Gin -TT -T " 



Sin Tra;^ 
sin (i 
sin TTX 

i cos 



sin 



sm -Kx 
sin 0(1 



cos0(^+a:o) 



sin Tra;^ 
cos 4){^-\-x\) 

1 C i n IT 'T" 1 



I sm Tra;^ 
i cos — ; 



(50) 

{x + {-lfx' + [-l)Hy)/2L. 

We will first reproduce the known results for the free 
Dirac fermions i.e. for clean graphene. To that end we 
rewrite the free energy (|4ip in the coordinate represen- 
tation as a function of 6: 



H<l>) 



w 



2 sin 0/2 



dxTr[V,Go(x,x,0)] . (51) 
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Substituting the bare Green function (1501) we obtain the 
generating function and the corresponding free energy 



W arcsin - 



■kL\/z 



(52) 



The corresponding distribution of transmission eigenval- 
ues then reads 



Po{T) = 



W 



1 



(53) 



The distribution P{T) is expected to be integrable and 
the integral gives the total number of open channels. 
However the integral of Eq. ([53]) diverges logarithmical 
at T = 0: we need to introduce a cutoff Tmin ~ e-^^^/"". 
However, the contribution of this cutoff to the higher 
moments is found to be exponentially small. Equa- 
tion (1531) coincides with the well-known result obtained 
by Dorokhov for the disordered metallic wiresi^ Thus, 
the transport of clean 2D Dirac fermions resembles the 
diffusive transport of non-relativistic electrons in quasi- 
one-dimensional systems in the presence of disorder. This 
nontrivial result can be explained by existence of the 
evanescent modes. 



C. The lowest order correction to free energy 

The perturbative corrections to the free energy of the 
system due to disorder can be expressed as a sum of loop 
diagrams without external legs. The lowest order con- 
tributions to the free energy are given by the one loop 
diagrams schematically shown in Fig. IHUa). The solid 
lines denote the propagator (|50l) in the presence of the 
boundary auxiliary field (f> and the dashed lines stands 
for disorder correlation functions. There are six topolog- 
ically equivalent diagrams with different disorder corre- 
lators corresponding to SR and LR correlated disorder 
of three types: random potential, random gauge (with 
two components x and y) and random mass. The corre- 
sponding integrals have the form 

^^SR,LR = 7^ / dWrV(r-r')Tr[E^G(r,r')S^G'(r',r)]. 



(54) 

with (?^(r), /i = 0,x,y,z given by Eq. (jH]) where we re- 
tain only the SR or LR part. These integrals diverge 
for r — > r', however, the divergent terms turn out to be 
(^-independent, and thus do not contribute to the phys- 
ical quantities (|47)) . Hence, it is convenient to consider 
the first derivative of the free energy with respect to (j) 
which is finite and determines the physical observables. 
The 0-dependent parts of diagrams with dashed line cor- 
responding to the three types of SR correlated disorder 
were computed in Ref . [2^ The result giving the linear in 
correction to the free energy of the clean sample ([5^ 
reads 



W(t> 

2^ 



("0 - Ctz) ■ 



(55) 




FIG. 4: (Color online) The one-loop LR disorder corrections 
to the generating function of transmissions (|62p for a — 1: 
ho{z) - solid black line and hz{z) - dashed green line. Inset. 
The one-loop LR disorder corrections to the free energy (|56p 
for a = 1: /o (<?!>) - solid black line and /z(0) - dashed green 
line. 



Note that the SR correlated random gauge potential 
does not contribute to the transport properties to one- 
loop order and this holds also to two-loop order. This 
is in agreement with the arguments of Ref. [13 that at 
zero energy the gauge potential can be eliminated by a 
pseudogauge transformation of the wave function. As a 
result, the transport properties are not infiuenced by ran- 
dom gauge potential despite the fact that it gives rise to a 
multifractal wave function \['(r) with a disorder strength 
dependent spectrum of multifractal exponentSi ^^i^'^" — 

The three diagrams with dashed line associated with 
correlation functions of the LR correlated disorder are 
computed in Appendix 1X1 The corresponding corrections 
to the free energy read 



W 



[fo{cl>)f3o - Mcf>)P, 



(56) 



We found that the LR correlated random gauge potential 
does not contribute to the transport properties to lowest 
order in disorder strength. The functions fp.{<f) are given 
by the following double integrals: 



dy / dc 







{ sinh y 


arctan ^ 







4 7r°-Ma;/sinh(y0/7r) 
(c2 + y^f^^ 

1— cose cosh y\ TT 

sin c sinh y J 2 



cosh y — cos c J ' 



(57) 



where the upper sign corresponds to fo{4') and the lower 
sign to fz{<f>)- The functions /^(i)f)) computed numerically 
for a particular value of a = 1 are shown in inset of Fig. |4| 
It is known that even in the case of uncorrelated disor- 
der the lowest order corrections to the density of states 
and the transport properties of 2D Dirac fermions are 
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2.0 FV 



1.5 




0.0 L 



0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

a 



FIG. 5: (Color online) The one-loop LR disorder corrections 
to the conductance (|64[) as functions of a: /o(0) - solid black 
line and /^(O) - dashed green line. 



insufficient. The leading corrections can be summed up 
with the help of the renormalization group methods that 
will be done in the next section. 



V. WEAK DISORDER RENORMALIZATION 
GROUP 

Straightforward dimensional analysis shows that the 
SR correlated disorder is dimensionless and hence 
marginally relevant in d — 2. The LR correlated dis- 
order is relevant in d = 2 for a < 2. In what follows it is 
convenient to introduce the rescaled disorder strengths: 
a{L) = a{L) and = /3(i)L'^-2. The lowest order 

corrections to the disorder strength and energy are given 
by the one-loop diagrams (b)-(e) shown in Fig. [31 The 
corresponding RG flow equations read 



dap 
dlnL 

d\nL 

da± 
d\nL 

dh 
ainL 

daz 
dlnL 

dXnL 

d\ae 
dlnL 



2<5o(ao + /3o + a_L + /3± + 5^ -t- /3^) 
(2 - a)/3o + 2;3o(5o + h + 



(58a) 

(58b) 
(58c) 



= 4(ao +/?o)(a^ +/?;.), 

= (2-a)^±, (58d) 

= -2^1 ~ 2a Jz + 2{&z + ao + /3o)(fi_L + 

-25^(ao+/3o), (58e) 
- (2 - a)/3, - 2/3^ - 2ci,/3, 



+2/3^(5_L+/3_L -5o-/3o), 



(58f) 




FIG. 6: (Color online) The LR correlated disorder correc- 
tions to the transmission eigenvalue distribution P{X) given 
by Eq. (pSj) for a — 1.5; po(A) - solid black line and Pz{X} - 
dashed green line as a function of A. Inset. The LR corre- 
lated disorder corrections to the Fano factor given by Eq. (|65p : 
/o"(0) - solid black line and /"'(O) - dashed green line, as func- 
tions of a. 



where we used the notation a± — ax + ay and — 
+ (3y Note, that in deriving the flow equations ([58]) 
we assume that 2 — a is small and perform 2 — a expansion 
similar to d—2 expansion in higher dimensions. In general 
in the presence of LR correlated disorder one has to rely 
on the double expansion in 2 — a and d — 2 similar to that 
for the 0^ model with correlated random bond disorder 
where one uses a double expansion in 4 — a and 4 — d at 
the upper critical dimension. 

The bare values of the disorder strengths and energy 
corresponding to the microscopic scale provide the initial 
condition for the RG equations (1551) . The renormalized 
disorder strengths, a{L), f3{L), and the energy e{L) ac- 
quire scale dependence on the ultraviolet cutoff length 
L. One has to stop the renormalization when either L 
reaches the system size or the energy e reaches the value 
of the cutoff A or the disorder strengths become of or- 
der one^SS. Once the renormalization has been done one 
can compute the observables by substituting the renor- 
malized quantities into the results of the perturbation 
theory. 

To renormalizc the corrections to the free energy ([551) 
and ([551) have to replace the bare coupling constants 
by the renormalized ones. As a result we obtain 



2ttL 
W 

2^ 



(59) 
(60) 



l + ao + fio + a±+P±+az+l3z, (58g) 



Thus, the SR correlated disorder does not modify the 
pseudodiffusive behavior to lowest order (fi ^ 0^) and 
the distributions of transmission eigenvalues is still given 
by the Dorokhov distribution (|53|) . The deviation from 
the pseudodiffusive regime can be found only to second 
order in disorder and the corresponding two-loop correc- 
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tions have the forro^i 

W(jP 



q21oop 



[(ao + Q!^)^wi((?!)) 



47rL 

+(ao + 3a2)(ao - az)]'^2(0), (61) 



where lui{4>) = const — ip((j)/Tr) — ip^—cjj/i:) and u)2{4>) = 
const + 7r^(0cot0 — ^)/(t>'^- Here ^"(2;) is the digamma 
function4^ 

On the contrary the LR correlated disorder leads to de- 
viation from pseudodiffusive transport already to lowest 
order in disorder. Indeed the renormalized corrections to 
the generating function and the distribution of transmis- 
sion eigenvalues read 



W 

W 
2^ 



(62) 
(63) 



Here hf^{z) = 2/^(2 arcsin-^i)/ sin </> and p^{X) = 
2Re/^(7r + 2iA)/7r. The functions h^{z) and p^(A) for 
particular values of a are shown in Figs. (jH) and ([5]). 
Note that the distribution can be used for direct cal- 
culation of transmissions moments for 1 < a < 2. This 
is different from the two-loop correction (|61l) due to SR 
correlated disorder found in Ref. the corresponding 
contributions to ■p(A) diverge at A = in a non-integrable 
way. This divergence has been attributed to the break- 
down of perturbative expansion in small disorder close to 
A = 0, i.e. for T k, 1. Nevertheless even for a < 1 one can 
compute the transport characteristics directly from the 
free energy. The correction to the conductance is given 
by 



G 



LR 



nhL 



/o(0)/3o-/^(0)/3. 



and the Fano factor can be written as 



F 



3 1 + ao - + /^/3o - nPz 



(64) 



(65) 



0=0 



/^(O) and //("(O) as functions of a are shown in Figs. ([5]) 
and (|6]), respectively. The case of generic disorder re- 
quires a numerical solution of the flow equations which 
strongly depends on the particular values of bare cou- 
plings. In the present paper we restrict our analysis to 
three cases when the system has only one type of disor- 
der: random scalar potential, random gauge potential or 
random mass disorder. 



A. Random scalar disorder 

Let us start the discussion of the different types of 
disorder with random scalar potential. In the presence of 
both SR and LR correlated scalar potential the solution 




FIG. 7: (Color online) The conductance as a function of a in 
the presence of random potential for L/lo — 0.5 - solid black 
line and in the presence of random mass - dashed green line. 



of the flow equations 



deep 
dlnL 

dlnL 
d\ne 
dlnL 



25o(ao + /3o), (66) 
(2-a)/3o + 2/3o(5o + /3o), (67) 
l + ao+/3n, (68) 



can be computed only numerically. However, we have 
found that the LR correlated disorder dominates over 
SR correlated disorder at large L for all bare disorder 
strengths such that a < (3. Since LR correlated disorder 
does not generate SR disorder itself we restrict ourselves 
to the case of pure LR disorder. Here and below we will 
measure the length in units of the bare ultraviolet cutoff 
given by vq/A. The solution of the flow equation ((67|) 
with the initial conditions i3o{l) = /3o = /3(A/z;o)''~^ 
reads 

/3o(i) = (a - 2) [2 + L--\a - 2(1 + • (69) 

The running disorder strength grows with the scale. At 
the Dirac cone the renormalization has to be stopped 
either at the scale of the system size or at the scale at 
which the disorder strength becomes of order unity. This 
scale computed from Eq. ([69)) reads 



2/3o 



(4 ~ a)/3o 



l/(2-a) 



(70) 



Iq is nothing but the zero-energy mean free path. For the 
system size L < Iq one can rewrite the running disorder 
strength in terms of Iq as follows 



ML) 



(2~a) 



{4-a){lo/L) 



2-a 



- 2 



(71) 



For finite energy the renormalization is limited by the 
scale at which the energy becomes of order of A. Sub- 
stitution of the solution (IMj) to the flow equation for the 
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energy (|68|) yields 

e{L) =eL[l- 2^o(i""' - 1)1 {a - 2)] 



-1/2 



(72) 



The renormalization stops when the running energy s{L) 
reaches the cutoff value A at the scale 



A 



L^{e) = 7 [1 - 2Mme\T-' - !)/(« - 2)] . (73) 

The competition between La and Zq introduces a new ex- 
ponentially small (in the limit a — > 2) in disorder energy 
scale r2-a given by equation iA(r2-a) = ^o, 




Ae 



+ 0{{2-af) 



.(74) 



For e S> T2-a the density of states can be found using 
the following scaling arguments. The running density of 
states approaches p = A/(27r-yg) at L = La- Taking into 
account that the density of states scales as pi = peL^ 
one can write the bare density of states as 



[l-2MA/k|)°-'-l)/(a-2)]-\(75) 



For £ < r2_a the density of states saturates at a finite 
value. This picture is in qualitative agreement with the 
prediction of SCBA computed in Sec. IIIII The results for 
SR correlated disorder case obtained in Ref. i23i can be 
reproduced by taking the limit of a — > 2. For instance, 
in this hmit we have Fq « A pl^^e~'^^'^^° . 

The conductance and the Fano factor in the ballistic 
regime L < Iq at the Dirac cone are given by 



G 

F 



1 + /^(O) ^0 (L) 



W 
nh L 
1 2 /^"(0)/3o(£) 
3 3l + /^(0)/3o(L)^ 



(76) 
(77) 



where fio{L) is given by Eq. (iTlT) . The conductance and 
the Fano factor computed for L/Iq = 0.5 are shown in 
Figs. JT]) and ^ as functions of a. The correction to 
the conductance due to LR correlated disorder in the 
ballistic regime is positive and increases with a while the 
correction to the Fano factor is small and negative. 



B. Random gauge potential 

We now turn to the case of random gauge potential. 
Inspired mostly by its relation to the Quantum Hall 
transitions,— this problem has previously motivated nu- 
merous studies of the multifractal spectrum for critical 
wave functions. -'^^^^^ Here we are mostly interested in the 
density of states and also transport properties of such 



Dirac fermions with correlated random gauge potential. 
In this case the flow equations reduce to 



9(5 _L 
aiuL 
Sine 







d\nL 



(2-a)/3i, 



(78) 
(79) 



The renormalized coupling constants have the trivial flow 
a±{L) = a± and I3±{L) — /Sj^L^"" where the bare dis- 
order strengths are a — a and I3± = /3±{A/vo)"'~'^- 
The LR disorder strength reaches unity at the scale 



lo = P 



-l/(2-a) 



Substituting the running couplings to 



the flow equation for the energy ([7^ we obtain 

£{L) = £Li+"^ exp \(3±{L^-'' - l)/(2 - a) 



(80) 



One has to stop renormalization at the scale La such 
that ^(La) = a. In the case of the system with only 
SR correlated random gauge disorder this scale is given 
by La^ = (e/A)~^^^, where we have introduced the dy- 
namic critical exponent z = 1 + a± . Note that this ex- 
ponent is non-universal and depends on the strength of 
disorder. In the presence of LR correlated disorder the 
cutoff scale computed up to subleading logarithmic cor- 
rections is 



-LR 



l-a, A 
-= — In — 

13± £ 



l/(2-a) 



(81) 



However, one has to stop renormalization at Zo for 
Iq < that introduces a new energy scale r2-a = 

Ae""'^/^^""-' which is exponentially small for a ^ 2. The 
bare density of states is then given by p = pe/ {eLj^) with 
p — A/ {27tvq ). Substituting the renormalized cutoff scale 
we obtain for the SR correlated disorder a non-universal 
power-law behavior: 



Psr(£) 



2ttv^ 



, z = l + aj_ (82) 



which was firstly derived in Ref. l27|. In the case of LR 
random gauge disorder we have 



PLR(e) 



A2 1 /2-a^^ Ay'/('""^ 



2ttvq e \ /3j_ £ j 

where in the last line we used the definition of the dimen- 
sionless disorder strength so that the dependence on the 
ultraviolet cutoff drops out from the density of states. 
Presumably Eq. (1551) is valid only for e > r2_a- In 
Sec. IVII we apply bosonisation technique to compute the 
density of states down to zero energy and show that the 
scaling behavior ((83)) actually holds up down to zero en- 
ergy. 

We have found above that the LR correlated random 
gauge disorder does not contribute to transport at the 
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Dirac cone. There are some general arguments that 
any random gauge potential cannot modify the trans- 
port properties. Let us first briefly recall the argument 
of Ref. [22!. To start, we consider the Hamiltonian (HJ-© 
with only a gauge field A: 



ivoa ■ V + * = 



(84) 

ft is known from vector analysis that any vector can be 
decomposed into the sum of a gradient and a rotational. 
Using that property, we can express the 2D vector A as 

l=Vx+(ixV)0. (85) 

Using this decomposition, we can rewrite a- Am the form 

a ■ A = a ■Vx + ^ ■ {z xV)(f>. (86) 

The mixed product a ■ {z x V0) = {a x z) ■ V0, and 
iaaz = [axz) so the Dirac equation ([5^ can be rewritten 
as: 

(87) 



Then, a pseudogauge transformation to a new wave func- 
tion ^ according to 



(88) 



turns the Dirac equation (|87p into the free Dirac equation 
without vector potential 



(89) 



and thus the transport properties of the initial model 1 
turn out to be the same as in the absence of the gauge 
potential. 

However, there are some subtleties in applying this ar- 
gument to correlated in space gauge potential. The dif- 
ficulty stems from the fact that the transformation ([55)) 
is not unitary. Indeed, if we denote the original and the 
transformed wave functions by 



v{x,y) 

Then, the normalization condition 



u{x,y) 
v{x,y) 



(90) 



1 = J dxdyMx,y)f + \v{x,y)\^], (91) 

transforms under pseudogauge transformation to 

1 = J dxdy[e-^'"''^''-y^\uix,y)f + e'^'"''^''-'y'^\vix,y)\'^]. 

(92) 

Therefore, the normalization condition (j92p is equivalent 
to the normalization condition (j9ip only for extended 
states, and only when 4'{x, y) is vanishing outside of a 




FIG. 8: (Color online) The Fano factor as a function of a in 
the presence of random potential for L/lo — 0.5 - solid black 
line (also Inset) and random mass - dashed green line. 



finite area. Indeed, in that case the normalization inte- 
gral is dominated by the asymptotic behavior of the ex- 
tended wavefunctions outside the finite area and is thus 
unchanged by the transformation. 

Thus, if (j>{x, y) is vanishing outside of a finite region, 
this would also imply a vanishing gradient and thus van- 
ishing correlations of the vector potential outside of this 
region. As a result, the correlations of the vector po- 
tential become necessarily finite ranged, in contradiction 
with the hypothesis of an infinite ranged power-law de- 
cay. This is in contrast with the case of a (5-correlated 
gauge potential which is compatible with a potential ex- 
isting only in a finite region of space. Nevertheless, we 
have not found any corrections to transport to one-loop 
order. 



C. Random mass disorder 

Let us now consider the system with only random mass 
disorder. The corresponding fiow equations are 



dlnL 
dp. 
dlnL 
dine 
dlnL 



-2&l^2aJ,, (93) 
(2 - a)/3, - 2/32 _ 2ci,/3„ (94) 
l + + (95) 



In the case of SR correlated disorder the running disorder 
strength approaches the Gaussian fixed point (a* = 0) 
so that the SR disorder is marginally irrelevant. This 
results in the logarithmic corrections to the scaling of 
the density of states 



Psr[£) = — 2"ln— • 



(96) 



In the most general case the fiow equations possess beside 
the unstable Gaussian fixed point (a* = /3* = 0) a non- 
trivial infrared stable fixed point (a* = 0, /3* = (2 — a)/2) 
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with eigenvalues Ai_2 = —2 + a negative for a < 2. The 
dynamic exponent describing the energy scahng is then 
given by 



z = 1 



/?: = 1 + (2 - a)/2. 



(97) 



The density of states has then universal scahng behavior 
Plr(£) - (98) 

The system of 2D Dirac fermions (or more precisely the 
pair of Majorana fermions) with random mass disorder 
is formally equivalent to two decoupled classical 2D Ising 
models with random bond disorder at criticality.^- It is 
known that uncorrelated random bond disorder is irrel- 
evant in RG sense resulting only in logarithmical cor- 
rections to the scaling of the pure Ising model. However, 
the LR correlated disorder is a relevant perturbation that 
changes the critical behavior?^ This latter result is in ac- 
cordance with our findings. 

The conductance and the Fano factor at the Dirac cone 
given by 



G 
F 



1 , 2 /f(0)(2-a)/2 



3 1-/^(0)(2-a)/2- 



(99) 
(100) 



The conductance and the Fano factor turn out to be 
also universal. They are shown in Figs. ^ and ([U as 
a functions of a. Since upon renormalization the dis- 
order couplings approach a fixed point of order 2 — a 
the system does not develop the mean free path scale. 
Thus, one can expect that the expressions for conduc- 
tance (|99l) and the Fano factor (jlOOp hold up to very 
large scale. Remarkably, that in contrast to uncorre- 
lated disorder which suppresses the Fano factor the cor- 
related disorder can enhance it. In the case of adatoms 
on the surface of topological insulator undergoing the 
paramagnetic-ferromagnetic Ising-like phase transition 
with rj — a — 1/4 one would expect on the basis of our 
one-loop treatment that the density of surface states be- 
haves as p{e) ~ e^/^^. Unfortunately, for a ^ 2 the low- 
est order correction in disorder becomes too large so that 
one cannot rely anymore on the one-loop approximation. 



VI. RANDOM GAUGE POTENTIAL: 
BOSONISATION 

In this section we reanalyze the problem of 2D Dirac 
fermions in the presence of LR correlated random gauge 
potential with the bosonisation technique. We will first 
give a detailed derivation of the bosonized action in the 
case of a general interaction, then discuss first the SR 
correlated disorder case^^i^M before turning to the LR 
correlated case and comparing our results with those of 
Sec. |Vl We start from the replicated action ([Ti)) with 
only terms with fi = x,y. In the partition function path 



integral, we introduce in the action the Matsubara time 
variable r = y/vo and we make the change of (indepen- 
dent) Grassmann variables according to; 



^a = ^. (101) 



The transformed action S = Sq + V , which we split for 
convenience into a free and interacting part, reads 

So = dxdr'^lpll-dr + iVoO^dx - Vn(Jy\i>a, (102) 
a 

V = 2^^o J dx dx dr dr' g{x — X ,T — t') 

\ a J \ b 

-(J2(^l^^)i^,r)] (Y.{yjUt){x',r')] ,(103) 



where we have dropped the tildes for clarity. The 
function g(x,T) for SR correlated disorder is given by 
g{x,T) — a±^5{x)6{vQT) and for LR correlated disorder 
by g{x,T) = /3±Aa{x'^ + Vqt'^)~'^/'^ . This action has the 
form of the action of a model of interacting fermions with 
an interaction that is non-local in Matsubara time.^^ The 
density of states can be calculated as 

p{e) = -- ImTrg(ii/„ ^ e + iO+) (104) 

TT 

with the trace of the Matsubara Green function given by 

TrGiii^n) = Tr[(it/„ + ^-i/)-i], 
d 

-[lnZ(zt'„)], 



d{iiyn) 



(105) 



where H is the corresponding Hamiltonian and Z the 
partition function. In Eq. (|105p , the average is taken with 
respect to the action (I102p - (|103p . It is convenient for 
the bosonization procedure to introduce the components: 



(106) 



and define: 

Jl=Y. Aa^L,a. Jr = Y1 Aa^R,a (107) 
o o 

to rewrite the interacting part of the action in the form: 



X [Jr(x, t)Jl{x\ t') + Jl{x, t)Jr{x',t')]. (108) 
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We can now apply the bosonization technique to the 
action (I102p - (|103p . First, the Hamihonian of the non- 
interacting part is rewritten in terms of the compo- 
nents (|106p as: 

Ho = J2^o,a (109) 

a 



(110) 



In bosonization, the fermion fields are expressed in terms 
of bosonic field o^^'^^ 9a and 0a as: 



1 



1 



(111) 
(112) 



with A a short distance cutoff, dxOa = Trlla and the fields 
(t)a and IIq satisfy the canonical commutation relations 
[(/)a(x), Iib{x')] — i5ab5{x — x'). The operators rjj^ji^ a are 
Majorana fermion operators that ensure anticommuta- 
tion of the fermion fields. The disorder-free part of the 
Hamiltonian has the bosonized form^ 



0,a 



dx 
"27rA 



dx cos 20a, 



(113) 



where we have chosen the same eigenvalue —i for all the 
products i]R,aVL,a- The disorder-free part of the action 
is then: 

So^iY^ J dxdTUadr(t)a - J drHo. (114) 



After integrating out the fields Ha in the path integral 
with action (jll4l) . the action of the sine-Gordon model is 
obtained,^'^ In the presence of disorder, we introduce the 
symmetric combinations of the bosonic fields: 



and the new fields 0a, II \ with 1 < A < ?i — 1 such that: 



n-l 



A=l 

ri-\ 



(116) 
(117) 



with: 



1 

n 



Ex A 

A=l 



6 — ^ab, 
n 



a=l 



(118) 



The conditions ()118p ensure that the new fields defined 
in Eq. (|116l) satisfy the canonical commutation relations. 
We can then express the disorder contribution (|108p to 
the action entierely in terms of He and 0c thanks to the 
relations: 



JH/L--^dxCl>c±^Uc. 

Zn Z 



(119) 



We will now discuss separately the two cases of SR and 
LR correlated disorder. 

SR correlated disorder. In the case of g{x,T) — 
a±5{x)S{voT), the disordered part of the action (|108p can 
be rewritten as 



V 



-2TTna±vo / dxdr 



(a.0c)' 



47r2 



(120) 



The fields He and IIa are then integrated out, leav- 
ing an action expressed purely in terms of 0c, 0a- The 
quadratic part of the action of the fields 0a is unchanged 
compared with the case without disorder, but the action 
of the field 0c becomes: 



dxdr 
2tt 



(5r0C 



vo{l - naj_) 



vq{1 + na^_){dx(t>cf 



.(121) 



The common scaling dimension of the fields cos 20a then 
becomes: 



dim. (cos 20a) = 1 



i^c-1 



where 



Kc 



1 — na^_ 
1 + ' 



and for n — >■ 0, 



dim. (cos 20a) — > 1 — a_L, 



(122) 



(123) 



(124) 



leading to the following renormalization group equation 

(125) 



dVn , X 



for Un. A strong coupling scale is reached for ~ 
\vnK/vo\-^'^'^+"^\ Using the Eq. pUSt and the scal- 
ing dimension (|124p the density of states is obtained 



as 



27.50.60 



Psr(£) 



1 



27rAi;o 



eA\ 1+° 



(126) 
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i.e a power-law enhancement with non- universal expo- 
nent is obtained. By comparing Eq. (|126p with the RG 
caculation result of Eq. (|82| we note that the two results 
are in perfect agreement provided the short distance cut- 
off is taken as A = wq/^- 

LR correlated disorder. In this case, introducing the 
Fourier transform g{q,uj) = ^il'^ + w^/wq)^''"^)/^ of 
g{x, r), we rewrite the action as: 



S = 



dqduj \(l)c{q,ujW 



27r2 
nvog 



2'KVr) 



1 - nvogiq,uj) 



, /" dxdr 

E y dxdT COS 2(/), 



vo 



+ {voqf (1+ 

(127) 



In general, a model with an action such as (|127p is not 
integrable. To estimate the free energy associated with 
(|127p we use the Gaussian Variational Method^i with 
(replica symmetric) variational action: 



dqduj \4)c{q,uj)\" 
27rwo 



27r2 



1 - nvog{q,uj) 



(voqf (1 + nvog{q,uj))] 



E 



2ttvo 



dxdr 
27r 



(9r(?!)A 



Wo 



1'o(9:e0a)^ 



and minimize the variational free energy: 

^var — Fq + {S — S'var)SvarJ 



Fn = -In 



(128) 

(129) 
(130) 



After some calculation, we find that: 



2 
Wo 



A 



1 



lim(02) = _ln 

n— >-0 z 



-nvo 



dujd{voq) UJ^ + {voq)"^ 



Kg, c.). (131) 



4^2 (w2 + („g^)2+^2)2- 

Solving the self-consistent equation (jl3ip . we obtain: 



C(a) V'^o 



a-2 



/3± 

C(«) 



1^ 

Wo 



a-2 



where Ty(x) is the already appeared in Sec. IIIII Lambert 
function**^ and we have introduced the function 



C(«) 



na{2 



— - sm — 

a) V 2 / 



(132) 



We find a density of states p{un ) 
for low energy as: 



I that behaves 



-2/(2-a) 



(133) 



hence, the density of state has a divergence for e — 
which is, however, integrable. Note that the re- 
sult (|133p is independent from the cutoff A. The density 
of states (I133P is expected to be valid down to zero en- 
ergy and agrees with the prediction of RG (|83| which is 
supposed to be valid at energies larger than the exponen- 
tially small in the limit a — >■ 2 energy scale. This proves 
that there exists only a single regime with the asymptotic 
behavior (|133p . 

Note, that the result (I133P differs from the density of 
states pie) ~ l/(eln|e|(6-'^)/(2-a)) obtained in Ref.M 
using a supersymmetric approach and a variational ap- 
proximation. We found that Eq. (20) of Ref. [H is not a 
correct solution of Eq. (19) of that paper. Upon finding 
the correct solution the subsequent calculations repro- 
duce our result (|133p . 



VII. CONCLUSIONS 

We have studied 2D Dirac fermions in the presence of 
LR correlated disorder with correlations decaying as a 
power law. In particular we have considered three types 
of disorder: random scalar potential, random gauge po- 
tential and random mass. Using the SCBA, weak dis- 
order RG and bosonisation technique we have computed 
the density of states modified by disorder in vicinity of 
the Dirac point of free fermions. Using a diagrammatic 
technique with matrix Green functions we have derived 
the full counting statistics of fermionic transport at low 
energy. Remarkably, in contrast with SR correlated dis- 
order the LR correlated disorder provides deviation from 
the pseudodiffusive transport already to lowest order in 
disorder. 

In the case of LR correlated random potential the pic- 
ture resembles that for SR correlated random potential. 
Using the SCBA and RG give qualitatively consistent 
picture: disorder generates an algebraically small energy 
scale below which the density of states saturates to a con- 
stant value while above this scale it is given by a corrected 
bare density of states. The correction to the conductance 
due to LR correlated disorder at the Dirac cone is posi- 
tive and increases with a while the correction to the Fano 
factor is small and negative. 

For LR correlated random gauge potential we have 
found that the density of states diverges at zero energy 
in an integrable way. This small energy behavior derived 
using bosonisation is completely consistent with the pre- 
diction of RG which is valid for larger energies. In partic- 
ular, the density of states is accessible in graphene using 
STM measurements that would allow one to measure the 
real exponent a describing the correlation of the random 
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gauge potential induced by ripples. We have found that 
the LR correlated random gauge potential does not con- 
tribute to the transport properties to one-loop order. 

In the case of LR correlated random mass disorder we 
have found a non-trivial infrared stable fixed point which 
controls the large scale properties of the disordered Dirac 
fermions. This results in a universal power law behavior 
of the density of states and universal transport proper- 
ties. Since the disorder couplings flow to the fixed point 
the system does not exhibit the mean free path scale. 
Thus, the conductivity and the Fano factor at the Dirac 
point are expected to have universal forms up to very 
large scales. Remarkably, that in contrast to uncorrelated 
disorder which suppresses the Fano factor the correlated 
random mass disorder enhances it. 
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Appendix A: One-loop diagrams contributing to the 
free energy 

In this appendix we compute the diagrams shown in 
Fig. EJa) with the dashed line corresponding to 3 differ- 
ent disorder correlators. To that end we substitute the 
bare Green function (ISUI) in Eq. (|5^ and evaluate the 



trace explicitly. Since the diagrams contain (j) indepen- 
dent divergent terms we will compute the derivatives of 
the diagrams with respect to (f). The diagrams with LR 
correlated scalar and random mass disorder then yield 



foA<f>) 



dy / dxi / dx 





X 



2Tr'^Aaysmh{y(f>) 



1 



cosh(7r?/) — cos(7r(a;i — 0:2)) 

I ), 

cosh(7ry) — cos(7r(a;i + 2:2)) / 



(Al) 



where the upper sign corresponds to /o and the lower sign 
to fz. The diagram with LR correlated random gauge 
disorder gives an expression which does not depend on 
(f) and thus it does not contribute to transport. We now 
change variables from xi and X2 such that cos(7r(a;i -|- 
X2)) ~ b and X2 — xi = c that formally can be written as 



1 1 

dxi / (ia;2/(cos[7r(a;i -I- a;2)], |x2 — xi|) 





1 cos TTC 



dc 



2db 



-J{b,c). 



, (A2) 

-1 

Applying transformation (|A2[) to Eq. (jAip and evaluating 
the integration over b we obtain Eq. ([57| . 
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